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Abstract. This paper proposes a dynamical notion of discrete geodesics, 
understood as straightest trajectories in discretized curved spacetime. 
The notion is generic, as it is formulated in terms of a general devia¬ 
tion function, but readily specializes to metric spaces such as discretized 
pseudo-riemannian manifolds. It is effective: an algorithm for computing 
these geodesics naturally follows, which allows numerical validation—as 
shown by computing the perihelion shift of a Mercury-like planet. It is 
consistent, in the continuum limit, with the standard notion of timelike 
geodesics in a pseudo-riemannian manifold. Whether the algorithm fits 
within the framework of cellular automata is discussed at length. 
Keywords: Discrete connection, parallel transport, general relativity, 
Regge calculus. 


1 Introduction 

Three reasonable hypotheses- bounded velocity of propagation of information, 
homogeneity in time and space, and bounded density of information -lead to 
the thesis that natural phenomena can be described and simulated by cellular 
automata. This implication has in fact been formalized into a theorem both in 
the classical J9] and the quantum case [T], albeit in flat space. Further evaluating 
this thesis leads to the project of selecting specific physical phenomena, such 
as gravitation, and attempting to describe them as cellular automata. A first 
step in this direction is to build discrete models of the phenomena. In the case 
of gravitation, this leads to the question we address in this paper: what is a 
discrete geodesics? 

Geodesics generalize the flat space notion of line, to curved spaces. A line is 
both the shortest, and the straightest path between two points, but in curved 
space the two criteria do not coincide [T2| . In computer graphics and discrete 
geometry, discrete geodesics as shortest path between two given point have been 
studied extensively mm- This is not the case of geodesics as straightest path 
given an initial point and velocity—with the noticeable exception of EE in the 
framework of simplicial complexes. Yet, it is this criterion that one must adopt 
in order to describe and simulate the timelike geodesics trajectories of particles. 
In this paper we adopt the dynamical, spacetime view on geodesics which is 
typical of numerical relativity EH- But instead of discretizing geodesics defined 




by partial differential equations in a continuous spacetime, we seek to discrete 
geodesics as a native notion of a discretized spacetime, for instance of a grid 
endowed with a metric. 

More precisely, the paper proposes both a notion of discrete-spacetime 
geodesics and a notion of discrete-time continuous-space geodesics (Section [2]). 
Both are generic, that is formulated in terms of a general deviation function, 
but readily specialize for metric spaces (Section [4]). They are effective: an al¬ 
gorithm for computing timelike geodesics naturally follows (Section [3]), which 
allows us to validate the notions numerically, by computing the perihelion shift 
of a Mercury-like planet (Section [5]). They are consistent with one another: the 
former is clearly a discretization the the latter. Moreover, the latter is proven 
to have the standard notion of continuous-spacetime geodesics in a Riemanian 
space as its limit, which validates both notions as legitimate discrete counter¬ 
parts (Section [9]). Whether the algorithm fits within the framework of cellular 
automata is discussed at length, as well as how this impacts on precision (Sec¬ 
tions [g|8]) . 

These results apply to natively discrete formulations of General Relativity 
such as Regge calculus ffZEl- For instance, our method yields perihelion shift 
computations of the right order, an issue in m pointed out in |3j. We discuss 
how our approach differs from m and why it fixes this issue. Finally, an often 
underestimated contribution is the pedagogical: the simple discrete model sum¬ 
marized in Figure [l] has continuum limit the complicated, well-known equations 
of ([9]) and ( |To| . 

Besides assessing this “digital physics” program, we believe that these results 
can be applied in any inherently discrete geometrical setting, in order to compute 
geodesics without the need to interpolate a continuous surface. Such applications 
may arise in computer vision and graphics |l4j including computer anatomy ma- 

2 Discrete geodesics 

Consider a discrete-time continuous-space spacetime Z x M n where Z is the 
discrete timeline and R n a continuous space. Consider a deviation function w 
from (Z x R ra ) 3 to R + , the number w(E, F, G) measuring how the path E, F, G 
deviates from “going straight ahead”. In this setting, a geodesic is a sequence 
of points in ZxR" {Ei)i such that for any i, w(£ l i _i, Ei, E i+ i) = 0. Such a 
property can be read as a condition on E l+ \: if the points 2£j_i and Ei are given, 
the geodesics must continue with a point E i+1 such that w(Ei_ 1; Ei, E i+ 1 ) = 0. 

Consider now a discrete spacetime M = Z x Z" where Z is the discrete 
timeline and Z" a discrete space and a deviation function w from M 3 to R + 
which, as above, measures how the path E , F , G deviates from going straight 
ahead. In this setting we cannot demand w(Ei_i, Ei, Ei+±) to be exactly zero, 
but we demand that it be a minimum with respect to spatial local variations of 
Ei+ 1 . 

Spatial local variations can be defined as follows. Let us write (xq,Xi, ... ,x n ), 
the coordinates of a point E in M, where xq is the time coordinate and x \,..., x n 



the space coordinates. Two points of M, (xq, x\,x n ) and (x' 0 , x ' x ,..., x' n ) are 
said to be spatial neighbors if Xq = x' 0 and for all i > 1, \x[ — Xi\ < 1. 

Thus, a discrete geodesics in M can be defined as a sequence of points in M, 
( Ei)i such that for any i, the deviation w(Ei-\,Ei,Ei + i) is a local minimum 
with respect to spatial local variations of E. L+ -\, that is for any spatial neighbor 
G of E i+ 1 , we have 


w{Ei-\, Ei, G) > w(Ei-i,Ei,E i+1 ) (1) 

Notice how this condition may be understood as a discrete counterpart of the 
Euler-Lagrange equation, in the spirit of El- 

3 An algorithm to compute a geodesic 

We now give a gradient descent-like algorithm to compute a discrete geodesic, 
(to, Ao), (ti, A\), (t 2 , A 2 ) given a deviation function w, a timeline to, ti, t 2 , •••, and 
two starting points A 0 an d A\. 

Assume, A,_j and A* are computed. To compute A, t+ i start with a point 
(U+i,C). Compute in((tj_i, Aj_i), (ti, A t ), (t i+ i, C')) for all 3” spatial neighbour 
C' of C. If they are all larger than w((ti_\, Aj_i), (ti, A), (t i+1 , C)) take C for 
Ai + 1 . Otherwise chose a C which minimizes w((tj_i, Ai_i), (ti, Aj), (tj+i, C')) 
and iterate, starting from this C'. 

Whether this iteration will eventually end depends, in general, on w (.,.,.). 
For instance, say that Aj_i), (ti, Aj), (t,+i, C')) increases as soon as 

AiC > ti +1 — ti. Then Aj+i will have to lie within distance £j +1 — £i of Aj, 
thereby imposing a bounded velocity c = 1, as well as enforcing termination. 


4 Distance induced deviation function 

Most of the times, the idea of deviating from going straight is induced from a 
notion of distance. Here is how. Suppose a distance function d, and define the 
three point distance function 

l(E, F, G) = d(E, F) + d(F, G). 

Intuitively, FG is understood to deviate from EF if it “leans” in some spatial 
direction FF', as witnessed by the fact that 

l(E,F',G) < l(E, F, G) 


for F' some neighbour of F. 

In a continuous-space discrete-time setting this would be formalized by letting 
w(E, (x 0 ,.. ,,x n ),G) be 


d 0 l(E, (x 0 , ..., x„), G )) 2 + ... + ( d n l(E , (xq, x n ),G)) 2 . 


( 2 ) 
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Fig. 1. Discrete geodesics seek to find G such that FG minimizes its deviation rel¬ 
ative to EF. In the case of metric spaces, FG is understood to “deviate towards FF' 
relative to EF”, whenever l(E, F 1 ,G) < l(E,F,G )—such deviations must be mini¬ 
mized. 

Thus, for continuous-space discrete-time geodesics (E^i each point Ei is a local 
extremum for Ei, E i+ i). 

In the discrete spacetime case, w(E, (xq, ..., x n ), G) is simply 
obtained by replacing partial derivatives with finite differences in Equation ([2|: 
d^l(E, (xq, ..., x n ), G ) becomes 



And, for discrete spacetime geodesics {Ei)i each point Ei minimizes the possibly 
non-zero w(Ei_ i, Ei, Ei + 1 ). 

5 Discrete Schwarzschild spacetime 

In this section, we give an example of discrete spacetime, which is a discretization 
of the Schwarzschild spacetime of General Relativity. 

Discretize spacetime down to A = 1 cm. Consider a star of mass M = 
2.10 3O /c<7 —alike the Sun. Its Schwarzschild radius is m = 2SM/c 2 = 3 km = 
3.10 5 cto. In order to evaluate distances, consider the metric tensor 



where r = \Jx 2 + y 2 , and let the distance function d be defined by 


d(E,F) = y / EF t s(£')EF. 


We study the geodesics trajectory of a planet, with respect to a timeline 0, r, 
2 t, 3r, ... with a = 10' and r = aA. Thus r = 10 'cm = 3.33.10 _4 s. The 
fake planet has parameters chosen so as to maximize relativistic effects: its first 










Fig. 2. The computed trajectory of a planet. 


point is E = {xq = 0,a;i = 10 8 cm = lOOOfcm,;^ = 0), and its initial velocity is 
vx = 0,vy = 2.10 -2 c = 6000/cto.s -1 . 

We compute the geodesics with respect to the w(., .,.) induced by d(., .) as in 
Section [4j and following the algorithm of Section [3j Recall that in this algorithm 
at iteration i the point A l+ \ is found by gradient descent starting from some 
point C. In the context of planetary movement, a good guess for C is obtained 
as follows. Define velocity S t = — Ai_ x and acceleration = S) — i, 

and make the guess that acceleration will remain constant, that is Ri+i = Ri. 
This would entail that A i+1 = Ai + Si + Ri, thus take C = A t + Si + Ri as the 
first guess and start exploring for the real Ai + 1 . Within reasonable ranges other 
heuristics- for instance, C = Ai + Si —lead to the same trajectories, but may 
require longer computation times. 

A run of the simulation is shown in Figures [2] and [3j Computation time is a 
few seconds. The code is available in [2]. The code is easily augmented to detect 
aphelion, typically a = lOOOfcTO, and perihelion, typically p = 150 km. 

The perihelion shift is visible on Figure [2] A well-known formula [5] states 
that perihelion shift in radians per revolution should be 

247T 3 L 2 67tS-M 37TTO 

a = T 2 c 2 (l - e 2 ) = c 2 L{ 1-e 2 ) = 

where T is the revolution period of the planet, L is the semi-major axis of the 
trajectory of the planet, e its eccentricity, and P = L( 1 — e 2 ) its parameter— 
recall that, by Kepler’s third law, T 2 = 47 t 2 L 3 /(SM), and m = 2SM/c 2 . Then 
an easy geometrical relation is 


1/a + 1 Ip 

hence 

a = (3/2)7rm(l/a + 1/p) 

which typically is 6.17 deg. The observed shift is around 6.27 deg. 









Perihelion 

t = 3580000000 cm x = -15031004 cm y = -1398397 cm 
angle = -174.6849818385271 deg 
distance = 1.5095913202506995E7 cm 
velocity = 0.13034045668302685 c 


Aphelion 

t = 7150000000 cm x = 99552205 cm y = 11035292 cm 

angle = 6.325378791430483 deg 

distance = 1.0016196478647615E8 cm 

velocity = 0.019968511905497616 c 

theoretical shift = 6.174380835177214 deg 

observed shift = 6.210787288401395 deg 


Perihelion 

t = 10730000000 cm x = -14847048 cm y = -2698767 cm 
angle = -169.69789787187233 deg 
distance = 1.5090333913952766E7 cm 
velocity = 0.13035029101632262 c 


Aphelion 

t = 14310000000 cm x = 97909925 cm y = 21872062 cm 

angle = 12.592542500595265 deg 

distance = 1.0032318032058926E8 cm 

velocity = 0.019937902185786747 c 

theoretical shift = 6.175065136027821 deg 

observed shift = 6.267163709164782 deg 


Fig. 3. Numerics of the computed trajectory of a planet. 


6 Cellular Automata in Mechanics 

As suggested in the Introduction, one motivation for discretizing General Rela¬ 
tivity is to describe the motion of a planet in a cellular automaton. 

Recall that, in the cellular automata vocabulary, a configuration cr is a func¬ 
tion which associates, to each cell C of the grid Z”, some internal state cr(C) 
taken in the set E. A cellular automaton is a function F from configurations to 
configurations, which has the following physics-inspired symmetries: 

— bounded velocity of propagation of information; 

— homogeneity time and space; 

— bounded density of information, that is E is finite. 

The state of a cell can be used to express the presence or the absence of a 
particle in this region of space. This way cellular automata can describe particle 
motions. For instance the simplest n-dimensional cellular automata—2 states, 






radius 1—can describe one particle motion among 3 n , as in each dimension, it 
could have velocity — 1, 0, or 1. 

To describe more complex motions, we must increase the number of states. 
For instance in a 1-dimensional automaton with radius 1, we can describe the 
motion of a particle that goes to the right at velocity 1/2, by alternating states 
si—stay still -and s 2 —step—, but also the motion of a particle that goes to the 
right at velocity 1, staying in the state S3. 

Another option is to increase the radius of the automaton. For instance, in a 
1-dimensional automaton, with radius 2 we can describe the motion of a particle 
that goes on the right at velocity 1 staying in a state Si or at velocity 2 staying 
in a state S2. Notice that modulo changing the units, the former behaviour can 
be obtained from the latter just by cell grouping. 

We want to address the following question: to what extent is the algorithm 
of Section [3] just a cellular automaton? 

7 Geodesics as Cellular Automata 

Whether the algorithm of Section [3] enforces a bounded velocity of propagation 
of information c = 1 depends, in general, on the properties of w (.,.,.). If such a 
velocity bound is enforced, then the motion of the body can be described in a 
cellular automaton of radius r. It is well-known that the velocity of a particle in 
a continuous Schwarzschild spacetime is bounded by c = 1. We conjecture that 
this is also the case for the discretized Schwarztchild spacetime. 

If w(., .,.) does not depend on space and time, then the algorithm clearly acts 
the same everywhere and everywhen, so that homogeneity is also enforced. In 
the important case where it depends upon a space-dependent metric, then this 
metric field has to be carried by the internal state of the cells, even if it does not 
contain a particle, so that homogeneity is still enforced. 

Let us evaluate whether bounded density of information holds. Even when 
w (.,.,.) does not depend on space and time, it is still the case that if a particle is 
at Ei, we need its velocity Ei_iEj to compute its next position Ei + \. But thanks 
to bounded velocity of propagation, and the fact that positions are discrete, the 
number of possible velocities is bounded above by 6 = (2a + l) n , so that bounded 
density of information is preserved. In the important case of a space-dependent 
metric carried by the internal state of the cells, whether bounded density holds 
depends upon whether we can assume that the metric field can be given with 
bounded precision. Even if this is not the case, notice that for a given cell, all that 
matters is to distinguish, for each input velocity of the particle, between b output 
candidate target cells. This map is a discrete counterpart to the connection 
associated to the metric. It contains just the finite amount of information that 
needs to be attached to the cell in order to compute geodesics. It could in fact 
be pre-compiled into each cell, thereby yielding a cellular automaton with 6 + 1 
internal states to code for presence and velocity, times b b to code for the discrete 
connection. 


8 Time versus space, precision 


Geodesics have been popularized by General Relativity. General Relativity likes 
to put space and time on an equal footing. Numerical schemes for General Rel¬ 
ativity ought to pursue that path, in particular it would be nice if the timeline 
of the computed geodesic were just 0, A, 2A, 3Z\, ... that is if a was equal to 1. 
In the scheme of Section [3j this choice leads to a cellular automaton of radius 
1, which is appealing, but it also restricts to b = 3 n the number of possible ve¬ 
locities. As we discussed in Section [6j this severely limits the number of motions 
that can be described. In the quantum setting, superpositions of basic velocities 
may compensate for this 033]. Classically, this is dramatic loss in precision. 
This is why, in Section [ 5 J we took a = 10 7 . 

However, we also saw that a radius of a' = 1 can be obtained from a cellular 
automaton of arbitrary radius a simply by grouping each a n hypercube of cells 
into one supercell. Each supercell now has an internal state in S' = S a . Notice 
that keeping the position of the single particle within the hypercube is crucial. 
Otherwise, all the velocities of norm less than one supercell are rounded up to 
the center of the supercell—and so the increased precision in the velocities is 
not much use. Hence, S' is really just coding for a velocity amongst b possibili¬ 
ties, which is appealing... but also for the position of the single particle within 
the hypercube, which perhaps is not so satisfactory. After all, what this space 
grouping has done is really just to hide the discrepancy between the discteti- 
zation step A and the computed geodesics timeline step a A, by hiding some of 
spatial precision within the internal space of the supercells. 

Hence, a » 1 appears to be fundamental requirement for precision. Notice 
that large values for a are better obtained by diminishing the discretization 
step A rather than augmenting the timeline step aA , as we cannot hope to 
achieve a pseudo-elliptic trajectory with just a handful of velocity changes per 
revolution. Running the simulations, it was indeed observed that a large, for 
instance a = 10', yields increased stability. But only to some extent: after a 
while the number of possible velocities b = (2 a + 1)" exceeds those which can 
be stored as a vector of machine-sized integers. 

It also helps to fine-grain the discretization step A, keeping a constant. Run¬ 
ning the simulations, it was indeed observed that this yields increased stability 
and convergence—at the expense of (reasonably) longer computation times. At 
some point, however, the finite-differences of <[3j) can become unstable, due to 
very small differences between l(E,F,G) and l(E,F',G) when FF' = 1, again 
hitting bounded machine floating point-arithmetic precision- but this can eas¬ 
ily be fixed by evaluating these derivatives with FF' a fraction of l(E,F,G) 
independent of A. 

9 Recovering continuous spacetime geodesics 

The algorithm of Section [3] is successful in computing geodesics in discrete time 
and space Z x Z n , in a way which is consistent with continuous-space discrete¬ 
time geodesics. We now explain how continuous-space discrete-time geodesics are 




themselves consistent with the standard geodesics of the fully continuous setting. 
For this question to make sense, we place ourselves in the case of Section [4j a 
distance-induced deviation function. 

As in Figure [4j consider three points E, F, G, taken at successive times 
0,r, 2 t. Let e be the distance EF, measured according to g(E), that is the 
proper time along EF and v = EF/e. In the same way, let e' be the distance 
FG according to g(F) and v' = FG/e'. 

We said that trajectory EFG is a continuous-space discrete-time geodesics 
if and only if it minimizes the distance EF + FG, with respect to infinitesimal 
changes of F into F'. Let us take FF' = Sd where d is a vector, normal with 
respect to g(F). 



Fig. 4. The continuum limit is obtained for e and 5 tending to zero. 


The distance EF 1 + F'G is given by: 

sj (ev + 5d)^g(E)(ev + Sd) + yJ (ev' — Sd)^ g{F'){ev' — Sd) 

Consider the first term. Its derivative with respect to 5 is 

d)g(E){ev + Sd) + (ev + 5d)t g(E)d 
2 \J (ev + Sdy g(E)(ev + Sd) 

Taken at S = 0 and using the symmetry of g(E) we get: 


<#g(E) (ev) 
y/(ev)ig(E)(ev) 


d ] g(E)v 


Consider the second term. If g(F') were just g(F), the same process would 
yield —d^g(F)v'. We would then just have <$g(E)v — $g(F)v’ = 0, yielding 
v' = g(F)~ 1 g(E)v. This is the equation derived in [T71 . and is in the same spirit 
as that obtained m in the framework of simplicial complexes. Unfortunately it 
does not yield accurate predictions for perihelion shift, as pointed out in H20 











and confirmed by our simulations. This is because one has to take into account 
that the variations of g(F) around F yield a third term: 


(eV)t(a g (F).rf)(eV) 

2 A /(eV) t 5(F)(eV) 


^ v'\dg{F).d)v' 


Let us emphasize that straightest geodesics on simplicial complexes [T5] do not 
see this term either: quite simply because a path EFG between two adjacent 
simplices sees the geometry of the first simplex—that is the term g(E )—and the 
geometry of the second simplex—that is the term g(F )—, but ignores the varia¬ 
tions of the geometry in some arbitrary direction FF'. In other words, simplices 
are usually thought of as polyhedrons of constant metric—but in order to be 
consistent with the continuum they must be interpreted as surfaces of constant 
metric derivatives. Altogether, we get that trajectory EFG is a geodesics if and 
only if 


d) (g(F)v' - g(E)v) = (dg{F).d)v' 


( 4 ) 


in every directions d. Unfortunately, this is still inconvenient to solve for v'. At 
this stage the traditional, continuous approach to geodesics follows two simpli¬ 
fying steps, which in the discrete setting translate into two approximations. The 
first step is to evaluate the condition for d only along the coordinate directions: 

(s(-F)aV - g(E)x.v) = | -v^{g(F) tX )v' (5) 

for all A. The second step is to realise that FG = EF + 0(e 2 ), and hence that 
v' = v + 0(e) and s' = e + 0(e 2 ). Indeed, in the continuum we assume that the 
time trajectory x{t) is twice diffentiable hence is has a first order expansion in 
t, from which we get FG = EF + 0(r 2 ) and EF ~ r^|, i.e. e ~ r. 

We get: 

g{F)\.v' = g{E) x.v + |v t (ff(F) jA )u. (6) 

V 1 = g(F)7x ( g{E)x-v + g(F),xvj ( 7 ) 

We get: 

gx■ v' = g X - v - egx-^v + | v^g iX v 
gx■ W - v) = | (u 1 #,xv - 2g Xv ^v ,l v v ) 
gx■ (v' - v) = | (g^xv^v 1 ' - gxv^v^v" - gx^^v") 

(v' -v) = -er ixv v ,l v v (8) 

where 

Ejj,v = g.x {gxv,[i T gx/j,,v g^u,x) /2 


( 9 ) 




Let us study the continuum limit e —> 0. Suppose E , F, G were three consecutive 
points along a curve x parametrized by its proper-time s. Since e was the proper 
time along EF, we have 


— = lim EF/e = lim v 
ds s —>-0 £—>-0 

and 

= lim(FG — EF)/e 2 = lim(u' — v)/e 
ds £— 

and hence 

d 2 x _ dx» dx v 

ds 2 ^ ds ds 

which is your traditional geodesics equation. 


( 10 ) 


10 Conclusion 

Summarizing, we have introduced a generic notion of discrete geodesics as straight - 
est trajectories in discretized spacetime, which are such that any three successive 
points E,F,G must minimize the deviation function w(E, F,G). Given E and 
F, G is implicitly determined: this can be viewed as a dynamical system and 
computed via a gradient descent algorithm. For a metric space, the canonical 
choice for w (.,.,.) measures how the length EF+FG varies with small variation 
of F. This was validated numerically, by computing the trajectory of a planet 
in discretized Schwarzschild spacetime, and recovering a perihelion shift of the 
right order. This was also validated by taking the continuum limit and recovering 
the standard geodesics equations on pseudo-riemannian manifolds. 

Part of our motivations were to evaluate the strength and limits of cellular 
automata. Recall that three well-accepted postulates about physics—bounded 
velocity of propagation of information, homogeneity in time and space, and 
bounded density of information- necessarily imply that physics may be cast 
in the framework of cellular automata—both in the classical and quantum set¬ 
tings m- Both theorems, however, rely on the implicit hypothesis of a flat 
spacetime. To which extent can cellular automata account for relativistic trajec¬ 
tories, that is geodesics? This paper shows that discrete geodesics can be cast in 
the framework of cellular automata, provided that a few extra assumptions are 
met: that the metric can be given with bounded precision, and that it has the 
property of fixing a velocity limit. These extra assumptions do not contradict 
the three postulates: they are but instances of them. 

Yet, this paper shows that a large discrepancy between the time discretiza¬ 
tion step the space discretization step is necessary in order to maintain a good 
precision on the velocity of particles. Namely, the number of particle velocities 
varies in (2 a + 1 ) n , with n the dimension of space and the radius a of the cel¬ 
lular automaton, which is therefore inherently large. Thus, computing discrete 
geodesics—and straight lines in euclidean space, for that matter—is local... but 


not that local. This may come as a surprise, and suggest that geodesics equations 
are better-behaved in the continuum. An alternative is to live with imprecise ve¬ 
locities. A planet is a collection of particles, and so it may be the average of their 
imprecise velocities which grants it a precise averaged velocity. In fact, a single 
particle is itself quantum, and may thus be in a superposition of these imprecise 
velocities, yielding a precise averaged velocity—as is made formal in the eikonal 
approximation [5]. This is in fact precisely what happens in quantum cellular 
automata models of quantum particles in curved spacetime, as shown in iazE|. 
All of these considerations suggest that nature’s way of working out timelike 
geodesics trajectories may in fact be emergent, from the simpler and more local 
behaviour of spinning particles. 

Hence, for future work, it may be interesting to look for discrete models 
based on spinning particles, oscillating along a few, cardinal, light-like direc¬ 
tions. These may in fact be closer to mimicking the real behaviour of fermions 
in curved spacetime, with the hope to recover the Mathisson-Papapetrou-Dixon 
equation—a generalization of the geodesics equation to spatially extended mas¬ 
sive spinning bodies—as emergent, in analogy with the continuum [5]. Such dis¬ 
crete models may be more local. Another approach is to work directly in terms 
of a discrete connection m- In the continuum, the Levi-Civita connection is 
axiomatized as being the unique metric-compatible and torsion-free connection. 
That given in Equation Q is exactly torsion-free, but interpreted as a discrete 
connection, as in Equation (|8|, it is metric-compatible only to first order. One 
can ask for both properties to be met exactly even in the discrete setting, this 
specifies the intersection of two ellipses. In 2-dimensions the number of solutions 
is finite, but this is not even the case in higher-dimensions: the axiomatization 
suggested by the continuum breaks down and demands fixing. 
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